Ground-state and dynamical properties of hard-core bosons on one-dimensional 
incommensurate optical lattices with harmonic trap 
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We study properties of the strongly repulsive Bose gas on one-dimensional incommensurate optical 
lattices with a harmonic trap, which can be deal with by using the exact numerical method through 
the Bose-Fermi mapping. We first exploit the phase transition of the hard-core bosons in the 
optical lattices from superfluid-to-Bose-glass phase as the strength of the incommensurate potential 
increases. Then we study the dynamical properties of the system after suddenly switching off the 
harmonic trap. We calculate the one-particle density matrices, momentum distributions, the natural 
orbitals and their occupations for both the static and dynamic systems. Our results indicate that 
the Bose-glass phase and superfluid phase display quite different properties and expansion dynamics. 
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I. INTRODUCTION 



Recently, the technical advances in cold atom trap- 
ping allowed the experimental realization of the Ander- 
son localization [l| in quantum matter waves \'& [3] . The 
good tunability and controllability of optical lattices of- 
fer myriad opportunities for studying the disorder effects 
in ultracold atom system So far, different techniques 
have been devised for the introduction of disorder in the 
ultracold atom system, such as speckle hied patterns in 
the optical lattice 0, H[ , random localized impurities by 
loading a mixture of two kinds of atoms with heavy and 
light masses @, and incommensurate bichromatic opti- 
cal lattices by superimposing two one-dimensional(lD) 
optical lattices with incommensurate frequency 0, 0] ■ I n 
particular, the experiment [7| has provided evidences of 
the existence of a Bose-glass (BG) phase [8|. 

Interesting phenomenons are expected to appear in dis- 
ordered systems when the interplay of disorder and in- 
teractions is taken into account. Interactions between 
atoms can be controllably tuned by Feshbach resonances 
in ultra-cold atom systems. While disorder can lead to 
localization of the wave function of a particle, delocal- 
ization can arise as the consequence of interactions in 
some many-body systems. Theoretically, for a repulsive 
Bose gas it has been predicted that there is a quantum 
phase transition from a superfluid phase to an insulating 
BG phase with localized single-particle states as disorder 
is increased , however unambiguous observation of 

the superfluid-Bose-glass transition is still under debate 
[H A lot of attention (l7M2l| has been payed to 

investigating the combined role of disorder and interac- 
tions in the strongly interacting ultra-cold atomic sys- 
tem. Except numerical or approximate approaches, the 
exact solution for the many-body systems with the in- 
terplay of disorder and interaction are rarely known. In 
this paper, we study the interacting bosons on the in- 
commensurate optical lattices with harmonic trap in the 
limit case with infinitely repulsive interaction which can 



be solved exactly. The ID Bose gas with infinitely repul- 
sive interaction is known as the hard-core boson (HCB) 
or Tonks- Girardeau (TG) gas (22[, which can be exactly 
solved via the Bose-Fermi map ping [22j and has attracted 
intensive theoretical attention [23- 25] . Experimental ac- 
cess to the required parameter regime has made the TG 
gas a physical reality [2(| [53] ■ For the HCB in ID optical 
lattices, it is convenient to use the exact numerical ap- 
proach proposed by Rigol and Muramatsu [28[. Follow- 
ing the exact numerical approach, we calculate the static 
properties of the hard-core bosons, such as one-particle 
density matrices, density profiles, momentum distribu- 
tions, natural orbitals and their occupations, to exploit 
the supcrfluid-to-BG phase transition for the systems in 
incommensurate optical lattices with harmonic confining 
trap. Furthermore, we study the nonequilibrium dynami- 
cal properties of expanding clouds of hard-core bosons on 
the ID incommensurate lattices after turning off the har- 
monic trap suddenly. We find that the expansion dynam- 
ics for the superfluid phase and BG phase exhibit quite 
different behaviors, which may serve as a signature for 
experimentally detecting the transition from supcrfluid- 
to-BG phase. 

The paper is organized as follows. In Section II, we 
present the model and the exact approach used in this 
paper. In Section III, we show properties of the ground- 
state for hard-core bosons on the incommensurate optical 
lattice with a harmonic trapping potential. Section IV is 
devoted to studying the nonequilibrium dynamics of the 
system after the harmonic trap is suddenly switched off. 
Finally, a summary is presented in Section V. 



II. MODEL AND METHOD 

In the present section we describe the exact approach 
which we used to study ID hard-core bosons on the in- 
commensurate lattice with an additional harmonic trap. 
Under the single-band tight-binding approximation, the 
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system of N hard-core bosons in the ID optical lattice 
can be described by the following Hamiltonian: 



h = -tj2( b l b i+i + H - c -) + E v ^ 



(i) 



where b\(pi) is the creation (annihilation) operator of the 
boson which fulfills the hard-core constraints [Hj], i.e., 
the on-site anticommutation {{bi,b\} — 1) and [£>»,&]] = 
for i ^ j; is the bosonic particle number operator; t 
is the hopping amplitude set to be the unit of the energy 
(t = 1); Vi is given by 



Vi = V/cos(a27ri + S) + V H {i - «o) S 



(2) 



Here Vj is the strength of incommensurate potential with 
a as an irrational number characterizing the degree of the 
incommensurability and 5 an arbitrary phase (in our cal- 
culation it is chosen to be zero for convenience, without 
loss of generality), Vh is the strength of harmonic trap 
and iq is the position of the vale of the harmonic trap. 

In order to calculate the properties of hard-core bosons, 
it is convenient to use the Jordan- Wigner transformation 
[29j (JWT) or Bose-Fermi mapping for the lattice model 



b \ = n e ~ inflh h =xi e+nflu u , (3) 



3=1 



13=1 



which maps the Hamiltonian of hard-core bosons into the 
Hamiltonian of noninteracting spinless fermions 

h f = - E(^+! + H - c -) + E V M> ( 4 ) 



where fj(fi) is the creation (annihilation) operator of the 
spinless fermion and n\ is the particle number operator. 
The ground-state wave function of the system with N 
spinless free fermions can be obtained by diagonalizing 
Eq.dU and can be represented as 



N 



i*£>=nE p *"tfi°>' 



(5) 



n— 1 i—l 



where L is the number of the lattice sites, N is the num- 
ber of fermions (same as bosons), and the coefficients 
Pi„ are the amplitude of the n-th single-particle eigen- 
function at the i-th site which can form an L x N matrix 

p m. 

In order to get the static properties of the ground- 
state, we calculate the one-particle Green function for 
the hard-core bosons defined by 



HCB 



bibles) = <* A |* B >, 



(6) 



where I^^^b) is the ground-state of hard-core bosons, 
and = (//n^e-^l^)) 1 , = 

/In^le'^^^l^i?)- Explicitly the state \^ A ) can 



be represented as 

Pi = -Pm for I < 
n<N,P A 
can get P B for the state 



n=l 2-iU 



= 1 and P A N+1 



1 P l if} |0) with 
= Pi n for I > i with 
= (I ^ i). Similarly we 
with the replacement of i 



1. 



pA 



by j. The Green function is a determinant dependent on 
the L x (N + 1) matrices P A and P B [H] 



G, 



(^ A \^ B ) = det (P A y P 



(7) 



The one-particle density matrix can be evaluated from 
the relation 



P,j 



(blb j } = G ij +5 ij (l-2Gu). 



(8) 



Alternatively one can use the method proposed by 
Pareder et al. [26] to calculate the one-particle density 
matrix directly. The momentum distribution is defined 
by the Fourier transform with respect to i — j of the one- 
particle density matrix with the form 



n(k) 



-ik(i-j) 



Pij, 



(9) 



where k denotes the momentum. The eigenfunctions or 
natural orbitals (4%) of the one-particle density matrix 
[30| can be obtained by solving 



L 



Pat] = \$ 



(10) 



which can be understood as being effective single particle 
states with occupations A^. For noninteracting bosons, 
all the particles occupy in the lowest natural orbital and 
bosons are in the BEC phase at zero temperature, how- 
ever only the quasi condensation exists for the ID hard- 
core bosons. For hard-core bosons in the incommensurate 
optical lattice, the strong incommensurate potential can 
destroy the quasi-BEC. 

The nonequilibrium dynamical properties of expanding 
clouds of hard-core bosons on ID incommensurate lattice 
after turning off the harmonic trap suddenly can be also 
calculated through the equal time Green function. With 
the JWT, we can express the equal time Green function 
for the expanding hard-core bosons as 



Gait) = ^ G HCB {t)\hb\\^ G HCB {t)) 



(ii) 



= (*f wi n j^'M} n e-**^\9$ (*)>, 

j3=\ 7=1 

where \^hcb(^)) * s the wave function of the hard-core 
bosons at t after turning off the harmonic trap and 
\^p(t)) is the corresponding one for the noninteracting 
fermions. The wave function \^p(t)} can be easily cal- 
culated with the initial wave function \^p) 
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FIG. 1: (Color online) The density profiles(a,c,e) and the local 
average density distributions(b,d,f) for the systems with 1000 
sites; 200 bosons; V H = 3 x 10~ 5 ; a = (a/5 - l)/2; and 
(a,b)Vr = 0, (c,d)Vi = 1, and (e,f)V> = 2.4. (g): The local 
average density distributions for the systems with 1000 sites, 
600 bosons, Vh = 3 X 10 -5 , a = - l)/2 and different Vi. 
(h): The mean particle fluctuation vs Vi for the systems with 
1000 sites, 200 bosons, Vh = 3 X 10~ 5 , a = (\/5- l)/2. 



which is still a product of time-dependent single-particle 
states, where h has been set to be unit of time in the 
evolution operator, H' F is the Hp in Eq.(j4|) with Vh = 0, 
and P(t) is the matrix of \$p(t)) in the same way as 
\%f). Then we can use the method described above to 
calculate the equal time Green function, yet we can get 
the one-particle density matrix, density profile, momen- 
tum distribution, natural orbitals and their occupations. 



III. STATIC PROPERTIES OF HARD-CORE 
BOSONS CONFINED IN HARMONIC TRAP 

Hard-core bosons in harmonic trap without the incom- 
mensuratepotential have been studied by Rigol and Mu- 
ramatsu [281 ] in detail. Quantitatively they characterize 
the system in a harmonic trap with the length scale set 
by combination lattice-confining potential £ = (Vn/t) 1 ^ 2 , 
and the associated characteristic density p = N/( [3lj |. 
They found that there is a critical characteristic density 
(p c ~ 2.6 — 2.7) in the system, for p < p c the whole sys- 
tem is in the superfluid phase at zero temperature, and 
for p > p c there is a phase separation in the system with 



a Mott insulating plateau in the middle of the trap with 
filling factor equaling one surrounding by the superfluid 
phase on the two sides. In this paper we mainly study 
the influence of the incommensurate potential which acts 
as the quasi-random potential and leads to localization 
of particles. We focus our study on systems with low 
characteristic density (jo < p c ). For system with high 
characteristic density the properties are similar to the 
one with low characteristic density except that there is a 
Mott insulator plateau in the system. The Mott insula- 
tor plateau is not dramatically influenced by the incom- 
mensurate potential because the particles in the Mott 
phase are already localized with one-particle density ma- 
trix pij — Sij . Without especially illustration, properties 
studied in this section are for systems with low charac- 
teristic density. 

Following the method described in the above section, 
we can get the one-particle density matrix. The den- 
sity profile is given by the diagonal elements of the one- 
particle density matrix, i.e., rii — pa. In the low charac- 
teristic density region, the density profiles for three differ- 
ent strength of the incommensurate potential are shown 
in FigJTJ In the absence of the incommensurate potential 
(Vi = 0), the density profile is shown in Fig. 1(a) with all 
the bosons being in the superfluid phase at the vale of the 
harmonic trap. As Vj increases but is still small, the den- 
sity profiles basically have the arc shape, but there are 
a lot of drastic oscillations in the profiles induced by the 
incommensurate potential. The amplitude of the oscil- 
lations becomes more and more large with the increase 
in Vi. When Vi becomes large enough, the density al- 
most oscillates in the range of to 1 . The density profile 
looks like a belt at the center of harmonic trap. Addi- 
tionally, the width of the density profile almost does not 
change with the change of Vj. For system in high char- 
acteristic density region, the Mott plateau in the density 
profile basically does not change with the increase in Vi 
because the particles in this phase are already localized. 
However the superfluid regions exhibit similar behavior 
as the system in low characteristic density region. 

Since there are many drastic oscillations in the density 
distributions, it is hard to tell that in which area the 
density is high in a density profile. In order to reduce the 
drastic oscillations, we define the local average density 
as 



M 

J2 n i+j /(2M- 

j=-M 



1), 



(13) 



where 2M + 1 is the length to count the local average 
density with M<L. In this work, we set M = 10. The 
local average density distributions for three different val- 
ues of Vi are also shown in FigQ] When Vj = 0, the 
local average density distribution is almost the same as 
the density profile with no drastic oscillations. As Vi in- 
creases, plateaus appear at the shoulders of the arc (see 
FigUJd)), then become wider and wider. Such kind of 
plateau is induced by the disorder, and thus we call it 
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Anderson plateau hereafter. When Vi > 2, the shape 
of the local average density does not change with the 
increase in Vj. In FigQ] (e) and (f), the density dis- 
tribution and the local average density distribution for 
Vi = 2.4 are displayed. The number of the plateaus and 
their locations are related to the particle number and the 
strength of the harmonic trap. In Fig|T] there happened 
to be only two plateaus. For system with high character- 
istic density, the local average density distributions for 
different Vj are shown in FigQJg). The Mott insulator 
plateau always exists in the trap center as Vj increases, 
except that the edges of the plateau are more and more 
dissolved into superfluid phase when Vj < 2. Meanwhile 
new plateaus appear at the two sides of the distribution, 
which is similar to the one with low characteristic density 
as Vi increases. The Mott plateau in the trap center is 
characterized by rii = 1 which does not oscillate against 
the particle density in the Anderson plateaus due to the 
existence of the incommensurate potential. 

The mean particle fluctuation for various value of Vi is 
shown in FigQ] (h). It is defined by n — 
Ki being the local particle fluctuation k.; = 



-yKi/L with 



iii — n'j for hard-core bosons. We can see that k de- 
creases as Vi increases because of the incommensurate 
lattice making the bosons difficult to hop. For system in 
high characteristic density region, the local particle fluc- 
tuation for particles in the Mott insulator plateau always 
equals zero as rt,; = 1 for any Vi, whereas the particle 
fluctuation is not equal to zero for particles in Anderson 
plateaus. So we can distinguish the Mott region from 
the superfluid region or Bose-glass region by the local 
particle fluctuation. 

In Figi2ta,b), we show the one-particle density matri- 
ces for different value of Vj and site. In the absence of 
the incommensurate potential (Vr = 0), the density ma- 
trix has a power-law decay with exponent —1/2 [2a ]. As 
Vi increases but is still small, the density matrices still 
have the power-law decay, but the exponents are smaller 
than —1/2 and have lots of oscillations induced by the 
incommensurate potential. Power-law decay of the den- 
sity matrix is the character of the system in superfluid 
phase, so when Vi is small, the system in the superfluid 
phase. But things change when the Anderson plateaus 
appear in the local average density distribution as Vj in- 
creases. The one-particle density matrices for the sites in 
Anderson plateaus have an exponential-law decay, while 
the ones for the sites out of the Anderson plateaus still 
have the power-law decay. Different from the case with- 
out the harmonic confining trap |l5j |. the- systems are 
not in a uniform phase. When Vi > 2, the local aver- 
age density distribution only consists of several Anderson 
plateaus, and the density matrix has the exponential-law 
decay which is the character of the system in Bose-glass 
phase. For system with p > p c , the one-particle density 
matrix is given by pij = Sij for i,j in the Mott insula- 
tor plateau, and the one for i,j on the shoulders in the 
density profile behaves like the density matrix of the sys- 
tem in low characteristic region. In Figj^tc) , we show the 
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FIG. 2: (Color online) One-particle density matrices for sys- 
tems with 1000 sites; 200 bosons; Vh = 3 X 10" 5 ; a = 
(V5- l)/2; and (a)i = 501, (b)i = 601. (c): The momen- 
tum distributions for systems with 1000 sites, 200 bosons, 
V H = 3 x 10" 5 , a={VE- l)/2. 



momentum distributions for different value of Vj . When 
Vi is small, the momentum distribution is similar to the 
case of Vi = and has a sharp peak at k = reflect- 
ing the coherence of the system. As Vi increases further, 
the peak becomes more and more shallow and the dis- 
tribution spreads wider. When Vr > 2, there is almost 
no obvious peak. The system is in the Bose-glass phase 
with all the effective single-particle states being localized 
states, and the coherence of the system becomes small. 

Now we study the properties of the natural orbitals and 
their occupations for the systems with the low character- 
istic density. The occupations of the natural orbitals for 
systems with different Vi are shown in FigO The occu- 
pations are plotted versus the orbital numbers w, and or- 
dered starting from the highest occupied one. When Vj is 
small, the occupation distribution has a sharp single peak 
at r\ — 1 which is the feature of bosons in contrast to the 
step function for fermions. With the increase in Vi, the 
occupation of the lowest natural orbital (Ai) decreases. 
When the system is in the Bose-glass phase for Vi > 2, 
no an obvious peak appears in the lowest natural orbital. 
Additionally we can find that there is a discontinuation 
in the occupation distribution at rj = N when Vj > 2. 
We use Z = \n — Xn+i [la] to characterize the discon- 
tinuation for a Bose gas with TV particles. In Figj3jd) Z 
versus Vi is plotted. There is an obvious change around 
Vi = 2. The difference of the occupations between sys- 
tems with the low and high characteristic density is that 
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FIG. 3: Occupations of the natural orbitals for systems with 
1000 lattice sites; 200 bosons; a = (a/5 - l)/2; Vh = 3 X 
10" 5 ; and (a)V = 0, (b)V = 1.2, and (c)V = 2.4. (d): The 
amplitude of the discontinuation (Z) which we define as Z = 
Ajv — Ajv+i relates to Vi for systems with 1000 lattice sites, 
200 bosons a = (a/5- l)/2 and V H = 3 X 10" 5 . (e,f): Profiles 
of the two lowest natural orbitals for systems with 1000 lattice 
sites, 200 bosons, a = (a/5 - l)/2 and Vh = 3 X 10" 5 . 



there is a plateau in the occupation for system in high p 
region due to the existence of a Mott insulator plateau. 
Such a plateau is basically not changed as Vi increases. 
In FigO we also show the effect of the incommensurate 
potential on the natural orbitals. We plot profiles of two 
lowest natural orbitals for different value of Vj. With- 
out the incommensurate lattice, the natural orbitals are 
similar to the ones of hard-core bosons in the harmonic 
trap without lattice. As Vj increases but is still small, 
the natural orbitals have similar shapes as the case of 
Vi = ones, but the widths of the wave functions be- 
come smaller and smaller due to the particles becoming 
more localized. When Vj > 2, all the natural orbitals 
spread only over a few lattice sites , an d they become 
typical Anderson localized state [HI, Hit. The system is 
in the Bose-glass phase with exponential-law decay one- 
particle density matrix when Vj > 2. For system with 
high p, the amplitude of the natural orbitals in the Mott 
insulator region has to be zero (see Ref. 28]). 



IV. DYNAMICAL PROPERTIES OF 
HARD-CORE BOSONS 

In this section, we study the nonequilibrium dynami- 
cal properties of expanding clouds of hard-core bosons on 
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FIG. 4: (Color online) The evolution of the density pro- 
file(a), momentum distribution(b), one-particle density ma- 
trix(c) with i = 501, occupations(d) of the natural orbitals for 
system with 1000 lattice sites, 100 bosons, a — (a/5 — l)/2, 
Vh — 3 x 10 -5 , and Vi = 0. t is the time after turning 
off the trap. The "fermions" curve is the momentum distri- 
bution for system with 1000 lattice sites, 100 free fermions, 
a = (a/5-1) A V H = 3xl0~ 5 , and Vi = 0. Insert of (d): The 
occupation of the lowest natural orbital vs t for system with 
1000 lattice sites, 100 bosons, a = (a/5- l)/2, Vh = 3 X 10" 5 , 
and Vi = 0. 



ID incommensurate lattice after turning off the harmonic 
trap suddenly. We first consider the situations that the 
systems arc in the supcrfiuid phase with small V/ before 
turning off the trap. The system with Vj = has been 
studied by Rigol and Muramatsu 



311 ] . In order to see 



clearly the effect of a nonzero value of V} and provide an 
example for comparison, we first present results for the 
case of Vi = following 31]. As shown in FigHJ the evo- 
lution of the density profile is ordinary, i.e., after turning 
off the harmonic trap the density profile spreads more 
and more wider as the time increases. As for the evo- 
lution of the momentum distribution, shortly after turn- 
ing off the trap, the peak at k = disappears. After a 
long time, the momentum distribution has similar shape 
as the distribution of the noninter acting fermions which 
does not change during the expansion 31]. Although the 
momentum distribution displays the behavior of fermion- 
ization |3lj . the modulus of one-particle density matrix 
still has the power-law decay which is the character of 
the bosons in superfluid phase. Here we note that the 
elements of the density matrix are complex numbers af- 
ter turning off trap. The occupations of natural orbitals 
basically remain unchanged, but the occupations of the 
lowest natural orbitals slightly increase during the ex- 
pansion. This can be understood as a result due to the 
increase in the size of the expanding atomic cloud. The 
evolution of the natural orbitals has similar behavior as 
the density profile which becomes more and more wider. 
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FIG. 5: (Color online) The evolution of the local average den- 
sity distribution(a), momentum distribution(b), one-particle 
density matrix with (c)i = 501, (d)i — 650, occupations(f) 
of the natural orbitals for system with 1000 lattice sites, 100 
bosons, a = (\Z5-l)/2, V H = 3x 10~ 5 , and Vi = 1. (e): The 
lowest natural orbitals at different time during expansion for 
systems with 1000 lattice sites, 100 bosons, a = (v5 — l)/2, 
Vb = 3 x 10" 5 , and Vi = 1. 



Now we consider the evolution of the system in the 
presence of the incommensurate potential. For the case 
that the system in a weak incommensurate potential is 
initially in the superfluid phase, the evolution of the den- 
sity profile is similar to the one with Vi = after turning 
off the harmonic trap. To give a concrete example, in 
Fig[5]we present the evolution of the local average den- 
sity for the system with Vi = 1. The density profile 
becomes more and more wider and the plateaus vanish 
during expansion. Comparing to the density profile in 
Fig 21 we can see that the expansion is much more slow 
as the incommensurate potential is added in which acts 
as a pseudo-random potential and prevents movement 
of particles. The evolution of the momentum distribu- 
tion is shown in FigjSJb). Shortly after turning off the 
harmonic trap, the peak at k — disappears. In the 
presence of the incommensurate potential, there appear 
some small peaks on the distribution. As time increases, 
the peaks of the momentum distribution become small, 
and the momentum distribution evolves to the one sim- 
ilar to the corresponding noninteracting fermions. The 
evolution of the one-particle density matrix is shown in 
FigJSfc,d). The density matrix \pij\ has the exponential- 
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FIG. 6: (Color online) The evolution of the local average den- 
sity distribution(a) and momentum distribution(b) for sys- 
tem with 1000 lattice sites, 100 bosons, a = (v§ — l)/2, 
V H = 3 X 10" 5 , and Vj = 2.5. 



law decay at t = for i in the Anderson plateaus. For a 
finite t, it still has the exponential-law decay despite the 
disappearance of Anderson plateau. For i out of the An- 
derson plateaus, \ptj\ has the power-law decay at t = 
, and it also has the exponential-law decay after turning 
off the harmonic trap. In FigEJf) we show the lowest 
natural orbitals for different time during expansion. The 
natural orbitals for system at different time have differ- 
ent shapes, which can be attributed to the scattering of 
initial wave function by the incommensurate potential. 
Correspondingly the occupations of the natural orbitals 
are different at different time during expansion. 

Finally, we study the evolution of system in the BG 
phase with Vj > 2. When Vj > 2, the system is in the 
BG phase with all the effective single-particle states be- 
ing Anderson-localized states. After turning off the har- 
monic trap, the evolution of the system is very slow, and 
basically remains unchanged even after a long time. As 
an example, we show the evolutions of the local average 
density and momentum distribution in Fig|51 in which 
we cannot find the difference between the initial one and 
the one after a long time. This clearly indicates that the 
system is completely pinned down by the strong incom- 
mensurate potential. The asymmetry of the local average 
density is caused by the asymmetry of the incommensu- 
rate lattice. 



V. SUMMARY 

In summary, we have studied the properties of hard- 
core bosons on incommensurate optical lattices with har- 
monic confining trap. Using the Bose-Fermi mapping and 
the exact numerical method proposed by Rigol and Mu- 
ramatsu (28l . [3lj . we calculate the one-particle density 
matrices, momentum distributions, natural orbitals and 
their occupations for both the static system and the dy- 
namic system at different time. Particularly, we exploit 
the phase transition from superfluid to the localized BG 
phase as the strength of the incommensurate potential 
increases from weak to strong, and the nonequilibrium 
dynamical properties of expanding clouds of hard-core 
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bosons on the incommensurate lattice after turning off 
the harmonic trap suddenly. For the optical lattice with 
harmonic confining trap, the density profiles show obvi- 
ous different characters for systems with weak, interme- 
diate and strong incommensurate potential. When the 
strength of the incommensurate potential becomes strong 
enough, Anderson plateaus are founded in the density 
distribution. The shape of the local average density pro- 
file shall not change with the increase in the strength of 
incommensurate potential if it exceeds a critical value of 
Vj = 2. When the harmonic trap is suddenly switched 
off, the expansion dynamics for the systems with Vj < 2 
and with Vj > 2 also exhibits quite different behaviors. 
All of these quantities give clear signature that there ex- 
ists a supcrfluid-to-Bose glass phase transition in the sys- 



tem when the strength of incommensurate potential ex- 
ceeds Vi = 2. Our study provides an exact static and dy- 
namic example which unambiguously exhibits the transi- 
tion from superfluid-to- Anderson-insulator in the incom- 
mensurate optical lattice. 
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